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I. INTRODUCTION 



proach has been successfully used in a time-dependent 
DP analysis of B° -> K + K~K° 0. 



The Cabibbo-Kobayashi-Maskawa (CKM) mecha- 
nism [l|, Q for quark mixing describes all transitions be- 
tween quarks in terms of only four parameters: three 
rotation angles and one irreducible phase. Consequently, 
the flavor sector of the Standard Model (SM) is highly 
predictive. One particularly interesting prediction is that 
mixing-induced CP asymmetries in decays governed by 
b —>■ qqs (q = u, d, s) transitions are, to a good approxi- 
mation, the same as those found in b — > ccs transitions. 
Since flavor changing neutral currents are forbidden at 
tree-level in the Standard Model, the b — ► s transition 
proceeds via loop diagrams (penguins), which are af- 
fected by new particles in many extensions of the SM. 

Various b — > s dominated charmless hadronic B de- 
cays have been studied in order to probe this predic- 
tion. The values of the mixing-induced CP asymmetry 
measured for each (quasi-)two-body mode can be com- 
pared to that measured in b — > ccs transitions (typi- 
cally using B° — > J/tpKg). A recent compilation [3[ 
of results shows that they tend to have central values 
below that for b — ► ccs. Recent theoretical evalua- 
tions 0, i, H 0, H H El [H d suggest that SM correc- 
tions to the b — > qqs mixing-induced CP violation param- 
eters should be small, in particular for the modes 4>K°, 
rj'K , and K S K S K S , and tend to increase the values, 
i.e. the opposite trend to that seen in the data. How- 
ever, there is currently no convincing evidence for new 
physics effects in these transitions. Clearly, more precise 
experimental results are required. 

The compilation given in Q includes several three- 
body modes, which may be used either by virtue of 
being CP eigenstates (K S K S K S , K s tt°tt ) [3 or be- 
cause their CP content can be determined experimen- 
tally (K+K-R ) [3, It also includes quasi-two- 
body (Q2B) modes, such as / (980)ijTg and p°(770)K s , 
which are reconstructed via their three-body final states 
(K s tt + it~ for these modes) . The precision of the Q2B ap- 
proach is limited as other structures in the phase space 
may cause interference with the resonances considered 
as signal. Therefore, more precise results can be ob- 
tained using a time-dependent amplitude analysis cov- 
ering the complete phase space, or Dalitz plot (DP), of 
B° — > iTg7r + 7r~. Furthermore the interference terms al- 
low the cosine of the effective weak phase difference in 
mixing and decay to be determined, helping to resolve 
ambiguities which arise from the Q2B analysis. This ap- 
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The discussion above assumes that the b — > s penguin 
amplitude dominates the decay. However, for each mode 
contributing to the K s ir + ir~ final state, there is also the 
possibility of a 6 — > u tree diagram. These are doubly 
CKM suppressed compared to the b — ► s penguin dia- 
gram (the tree is 0(A 4 ) whereas the penguin is Q(\ 2 ), 
where A is the usual Wolfenstein parameter [H [TtJ )• 
However, hadronic factors may enhance the tree ampli- 
tudes, resulting in a significant "tree pollution." These 
hadronic factors may be different for each Q2B state, 
thus the relative magnitudes of each tree and penguin 
amplitudes, \T/P\, and the strong phase difference may 
be different as well. Nontheless, the relative weak phase 
between these two amplitudes is the same - and in the 
Standard Model is equal to the CKM unitarity triangle 
angle 7. An amplitude analysis, in contrast to a Q2B 
analysis, yields sufficient information to extract relative 
phases and magnitudes. Measurements of decay ampli- 
tudes in the DP analysis of B° — > K s ir + ir~ (and similar 
modes) can therefore be used to set constraints on the 
CKM parameters (p, rf) [H, [H M, HI • 



Recently published results on time-dependent DP anal- 
ysis of B a — > K s tt + it~ are available 22j. Previous 



studies of the B° — ► KqTt + tt~ decay were either based 
on a Q2B approach [23|, or were amplitude analyses 
that did not take into account either time-dependence 
or flavor-tag dependence [24j |. The available results for 
B° — > Kgir + ir~ are consistent with studies obtained from 
other B — > Kirir decay modes: K + ir~Tr° [25|, |26| and 
K + n + n~ [23, [28|. The latter results indicate evidence 
for direct CP violation in the B + — > p°(770)K + channel. 
If confirmed, this will be the first observation of CP viola- 
tion in the decay of any charged particle. The relevance 
of B — ► Kitty is further highlighted by recent theoreti- 
cal calculations [2!| suggesting that large CP violation 
effects are expected in several B — > K*ir and B — > Kp 
resonant modes. 



In this paper we present results from a time-dependent 
amplitude analysis of the B Q — > Kgir + ir~ decay. In 
Sec. |H] we describe the time-dependent DP formalism, 
and introduce the signal parameters that are extracted 
in the fit to data. In Sec. lIIII we briefly describe the BaBAR 
detector and the data set. In Sec. II VI we explain the se- 
lection requirements used to obtain the signal candidates 
and suppress backgrounds. In Sec. [V] we describe the fit 
method and the approach used to control experimental 
effects such as resolution. In Sec. [VI] we present the re- 
sults of the fit, and extract parameters relevant to the 
contributing intermediate resonant states. In Sec. IVIII 
we discuss systematic uncertainties in the results, and 
finally we summarize the results in Sec. IVIIII 
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II. ANALYSIS OVERVIEW 

Taking advantage of the interference pattern in the DP, 
we measure relative magnitudes and phases for the dif- 
ferent resonant decay modes using a maximum-likelihood 
fit. Below, we detail the formalism used in the present 
analysis. 



• to is the invariant mass of the decay products of 
the resonance, 

• Rj(m) is the resonance mass term or "lineshape" 
(e.g. Breit-Wigner) , 

• L is the orbital angular momentum between the 
resonance and the bachelor particle, 



A. Decay amplitudes 

We consider the decay of a spin-zero B° with four- 
momentum pb into the three daughters tt + , 7r~, and K s 
with p +1 p-, and po their corresponding four-momenta. 
Using as independent (Mandelstam) variables the invari- 
ant squared masses 



m K° s -K+ 



= fa 
= (P- 



Po) 



(1) 



the invariant squared mass so = m 2 +7T _ = (p+ + P-) 2 
can be obtained from energy and momentum conserva- 
tion: 



so 



m B<> + 2m 2 r+ + m 2 K Q 



s+ - a- ■ (2) 

The differential B° decay width with respect to the vari- 
ables defined in Eq. ([T|) (i.e. the Dalitz plot) reads 

\A\ 2 



1 



ds+ds- 



(3) 



(2tt) 3 32m 3 BO 

where A is the Lorentz-invariant amplitude of the three- 
body decay. In the following, the amplitudes A and 
A correspond to the transitions B° — > K s ir + Tr~ and 
B° — > K s ir + ir~, respectively. We describe the distribu- 
tion of signal events in the DP using an isobar approx- 
imation, which models the total amplitude as resulting 
from a coherent sum of amplitudes from the N individual 
decay channels 



A(s A 



A(s^ 



N 
N 

E 

3=1 



CjFj{s+,S- 



(4) 



(5) 



where Fj are DP-dependent dynamical amplitudes de- 
scribed below, and Cj complex coefficients describing the 
relative magnitude and phase of the different decay chan- 
nels. All the weak phase dependence is contained in Cj, 
and Fj contains strong dynamics only; therefore, 

F.,- («+,«_) - F,( s _, s+ ) . (6) 

The resonance dynamics are contained within the Fj 
terms, which are represented by the product of the in- 
variant mass and angular distribution probabilities, i.e., 



Rj(m)X L (\p*\ r')X L (\q\ r)Tj(L,p, q) 



(7) 



where 



• p* is the momentum of the bachelor particle eval- 
uated in the rest frame of the B, 

• p and q are the momenta of the bachelor particle 
and one of the resonance daughters, respectively, 
both evaluated in the rest frame of the resonance 
(for KgTT~, Kgir + , and tt + tt~ resonances, q is as- 
signed to the momentum of the Kg, ir + , and 7r~, 
respectively) , 

• Xl are Blatt-Weisskopf barrier factors [3fJ with 
parameters r' (taken to be 2(GeV/c)~ 1 ) and r 
(given in Table [TJ> , and 



• Tj{L,p,q) is the angular distribution: 



L = : Tj = 1 , 

L = 1 : Tj=-4p-q, 



(8) 
(9) 



L = 2 : T 3 = -[i{p-qf-{\p\\q\f] ■ (10) 



The helicity angle of a resonance is defined as the an- 
gle between p and q. Explicitly, for Kg%~, K s ir + , and 
7r + 7r~ resonances the helicity angle is defined between 
the momenta of the bachelor particle and of the K s , 7r + , 
and 7T~ , respectively, in the resonance rest frame. 

For most resonances in this analysis the Rj are taken 
to be relativistic Breit-Wigner (RBW) [3l| lineshapes: 



Rj(m) 



(ttt-q — to 2 ) — imoT(m) ' 



(11) 



where mo is the nominal mass of the resonance and T(m) 
is the mass-dependent width. In the general case of a 
spin- J resonance, the latter can be expressed as 



r(m) = To 



/ \ 2J+1 










V TO / 



(12) 



The symbol Tq denotes the nominal width of the reso- 
nance. The values of too and Tq are listed in Table [J 
The symbol qo denotes the value of q when to = Too- 

For the /o(980) lineshape the Flatte form [32j is used. 
In this case the mass-dependent width is given by the 
sum of the widths in the mr and KK systems: 



T(to) 



T^{m) + T KK (m), 



(13) 
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where 



IVn-(m) = g 



(14) 



-^/l - 4m^ ± /m 2 



r/f/f(m) = 5if ( -4m|. ± /m 2 + 



(15) 



iyd - 4m| r0 /m 2 



The fractional coefficients arise from isospin conservation 
and g n and qk are coupling constants for which the values 
are given in Table [T] 

For the p° (770) we use the Gounaris-Sakurai (GS) pa- 
rameterization [33| , that describes the P-wave scattering 
amplitude for a broad resonance, decaying to two pions: 



where 

f(m) = T 



1 + d ■ T /m 



(to 2 , — to 2 ) + /(to) — i moT(m) 



Q 3 o 



q ( h(m) — h(m Q ) ) + 
dh 



(16) 



(17) 



TO, 







TO 



So 2 



dm 2 



and the function h(m) is defined as 



with 
dh 



dm 2 



him) — — — In 

7T TO 



2g 



2TO7T 



2m 2 , 



27TTO5 



(18) 



(19) 



The normalization condition at Rj (0) fixes the parameter 
d = /(0)/(roTOo). It is found to be: 



d - - ™k In ( m ° + 2q ° 

V 2m v 



7T <?o 



27rg 



mlm 



(20) 



The + component of the if 7r spectrum is not well un- 
derstood [H, HH; we dub this component (Kir)^ and 
use the LASS parameterization [34[ which consists of the 
if* (1430) resonance together with an effective range non- 
resonant (NR) component: 



Rj(m) 



mKrr 



gcot ( 

„2iS B 



iq 



(21) 



(mg - m 2 K7T ) - im T 



mKTT qo 



where cot Sb — ^ + \rq. The values we have used for 
the scattering length (a) and effective range (r) param- 
eters of this distribution are given in Table Q] The effec- 
tive range part of the amplitude is cut off at m™£ off = 



1800 MeV/c 2 . Integrating separately the resonant part, 
the effective range part, and the coherent sum we find 
that the iV*(1430) resonance accounts for 81.7%, the ef- 
fective range term 44.1%, and destructive interference be- 
tween the two terms is responsible for the excess 25.8%. 

A flat phase space term has been included in the signal 
model to account for NR B° — > Kgir + ir~ decays. 

We determine a nominal signal Dalitz^plot model us- 
ing information from previous studies [23l [24| and the 
change in the fit likelihood value observed when omitting 
or adding resonances. The components of the nominal 
signal model are summarized in Table [U Other compo- 
nents, taken into account only to estimate the DP model 
uncertainty, are discussed in Sec. I VIII 



TABLE I: Parameters of the DP model used in the fit. Values 
are given in MeV(/c 2 ), unless mentioned otherwise. The mass 
and width for the fx (1300) are averaged from results in B + — » 
K+TV--K+ Dalitz analyses [27|, [H. 



Resonance 


Parameters 


Lineshape 


Ref. for 

Parameters 


/o(980) 


mo = 965 ± 10 
g n = 165 ± 18 
g K = 695 ± 93 


Flatte 


[36] 


p°(770) 


m = 775.5 ± 0.4 
T = 146.4 ± 1.1 

r = 5.3«^(GeV/c)- 1 


GS 


121] 


iT" + (892) 
K— (892) 


m = 891.66 ±0.26 
To = 50.8 ± 0.9 
r = 3.6 ± 0.6 ( GeV/cT 1 


RBW 


[31] 




m = 1415 ± 3 
r = 300 ± 6 

m^ off = 1800 
a = 2.07±0.10(GeV/c)- 1 
r = 3.32 ± 0.34 (GeV/c)- 1 


LASS 


[27] 


/ 2 (1270) 


rn = 1275.4 ± 1.1 

To = 185.2t^ 
r = 3.0(GeV/c)- 1 


RBW 


[3JJ 


/x(1300) 


mo = 1471 ± 7 
To = 97 ± 15 


RBW 


[27, 28] 


NR decays 




flat phase space 




Xrf) 


mo = 3414.75 ± 0.35 
To = 10.4 ±0.7 


RBW 


[31] 



B. Time dependence 

With At = i s i g — t tag defined as the proper time 
interval between the decay of the fully reconstructed 
B° — > Kg7r + ir~ (^si g ) and that of the other meson 
(£?" ag ) from the T(AS), the time-dependent decay rate 
\A+ e (At)\ 2 (\A; ig (At)\ 2 ) when the B t ° ag is a 5° (S°) is 
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^iven by 



IAl g (A*)| 2 



|At|/r B o 



4r B o 

|2 



|^| 2 + |^| 2 

T (|^| 2 - \A\ 2 ) cos(Aro d At) 



± 21m [AA*] sin(Am d At) 



(22) 



where t^o is the neutral B meson lifetime and Ato^ is the 
B°B° mass difference. In the last formula and in the fol- 
lowing, the DP dependence of the amplitudes is implicit. 
Here, we have assumed that there is no CP violation in 
mixing, and have used a convention whereby the phase 
from B°B° mixing is absorbed into the B° decay ampli- 
tude (i.e. into the Cj terms). In other words, we assume 
that the B°B° mixing parameters satisfy \q/p\ = 1 and 
absorb q/p into Cj. Lifetime differences in the neutral B 
meson system are assumed to be negligible. 



storage ring at SLAC between October 1999 and August 
2006. The sample consists of an integrated luminosity of 
347.3 fb~\ corresponding to (383 ± 3) x 10 6 BB pairs 
collected at the T(4S) resonance ( "on-resonance" ) , and 
36.6 fb -1 collected about 40 MeV below the T(4S*) ("off- 
resonance" ) . 

A detailed description of the BABAR detector is pre- 
sented in Ref. [38|. The tracking system used for track 
and vertex reconstruction has two components: a silicon 
vertex tracker (SVT) and a drift chamber (DCH), both 
operating within a 1.5 T magnetic field generated by a 
superconducting solenoidal magnet. Photons are iden- 
tified in an electromagnetic calorimeter (EMC). It sur- 
rounds a detector of internally reflected Cherenkov light 
(DIRC), which associates Cherenkov photons with tracks 
for particle identification. Muon candidates are identified 
with the use of the instrumented flux return (IFR) of the 
solenoid. 



C. The square Dalitz plot 



IV. EVENT SELECTION AND BACKGROUNDS 



Both the signal events and the combinatorial e + e~ — > 
qq (q = u, d, s, c) continuum background events populate 
the kinematic boundaries of the DP due to the low final 
state masses compared with the B° mass. The represen- 
tation in Eq. ([3]) is inconvenient when empirical reference 
shapes are to be used. Large variations occurring in small 
areas of the DP are very difficult to describe in detail. We 
therefore apply the transformation 



dsj- ds- 



| det J\ dm' d9' 



(23) 



which defines the square Dalitz plot (SDP). The new co- 
ordinates are 



1 



in 



■ arccos 



m 



-. max 
l 



- 1 



1 



(24) 



where mo 



m B o 



so is the 7T + 7r 
and m™ m = 



invariant mass, m™ &yi 
■ the kinematic limits 



of $o is the tt + tt~ resonance helicity angle and J 
is the Jacobian of the transformation. Both variables 
range between and 1 . The determinant of the Jacobian 
is given by 

i , ti A | * I, *< 9m dcosd . . 

|detJ| = 4|p + ||p |m -^ 7 -^ r - , (25) 



89' 



E*^-ml + and |p5| = 



E* 2 



where | P ;| = „ v+ _ lt . 0l - ...^ 

and where the n + (Kg) energy E5_ (Eq), is defined in the 
7r + 7r~ rest frame. This transformation was introduced 
in Ref. [3?|], and has been used in several B decay DP 
analyses. 



III. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 



We reconstruct B° — > KgTr + ir~ candidates from pairs 
of oppositely-charged tracks and a Kg — > tt + it~ candi- 
date, which are required to form a good quality vertex. In 
order to ensure that all events are within the DP bound- 
aries, we constrain the invariant mass of the final state 
to the B mass. For the ir + ir~ pair from the B, we use 
information from the tracking system, EMC, and DIRC 
to remove tracks consistent with electron, kaon, and pro- 
ton hypotheses. In addition we require at least one track 
to be inconsistent with the muon hypothesis based on 
information from the IFR. The Kg candidate is required 
to have a mass within 15 MeV/ c 2 of the nominal K° 
mass [3l| , and a lifetime significance of at least five stan- 
dard deviations. The last requirement ensures that the 
decay vertices of the B° and the Kg are well separated. 
In addition, combinatorial background is suppressed by 
requiring the cosine of the angle between the Kg flight 
direction and the vector connecting the B-daughter pions 
and the Kg vertices to be greater than 0.999. 

A B-meson candidate is characterized kinemat- 
ically by the energy-substituted mass Toes = 



y/( S /2 + Pi -p B )yE? 



Pi 



and energy difference 



AE = E* B — \\fs~, where (E b ,Pb) and (J? i( pj) are 
the four-vectors of the B-candidate and the initial 
electron-positron system, respectively. The asterisk 
denotes the T(AS) frame, and s is the square of the in- 
variant mass of the electron-positron system. We require 
5.272 < m ES < 5.286 GeV/c 2 and \AE\ < 0.065 GeV. 
Following the calculation of these kinematic variables, 
each of the B candidates is refitted with its mass 
constrained to the world average value of the B-meson 
mass (3l| in order to improve the DP position resolution, 
and ensure that Eq. |(2J) holds. 

Backgrounds arise primarily from random combina- 
tions in continuum events. To enhance discrimination 
between signal and continuum, we use a neural network 
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(NN) [39| to combine four discriminating variables: the 
angles with respect to the beam axis of the B momen- 
tum and B thrust axis in the T{AS) frame, and the ze- 
roth and second order monomials L0.2 of the energy flow 
about the B thrust axis. The monomials are defined by 
L n = J2iPi x |cos#j| , where 6i is the angle with respect 
to the B thrust axis of track or neutral cluster i and pi 
is the magnitude of its momentum. The sum excludes 
the B candidate and all quantities are calculated in the 
T(4S) frame. The NN is trained using off-resonance data 
as well as simulated signal events, all of which passed the 
selection criteria. The final sample of signal candidates 
is selected with a requirement on the NN output that re- 
tains 90% of the signal and rejects 71% of the continuum. 

The time difference At is obtained from the measured 
distance between the positions of the B® ig and S° ag decay 
vertices, using the boost f3j = 0.56 of the e + e~ system. 
B° candidates with |Ai| > 20 ps are rejected, as are can- 
didates for which the error on At is higher than 2.5 ps. To 
determine the flavor of -B° ag we use the B flavor tagging 
algorithm of Ref. [4(j. This algorithm combines several 
different signatures, such as charges, momenta, and decay 
angles of charged particles in the event to achieve optimal 
separation between the two B flavors. This produces six 
mutually exclusive tagging categories: lepton, two dif- 
ferent kaon categories, slow pion, kaon-slow pion, and a 
category that uses a combination of other signatures. We 
also retain untagged events in a seventh category since 
although these events do not contribute to the measure- 
ment of the time-dependent CP asymmetries they do pro- 
vide additional statistics for the measurements of direct 
CP violation and CP-conserving quantities such as the 
branching fractions 4l|. Multiple B candidates passing 
the full selection occur between ~ 1% of the time for NR 
signal events and ~ 8% of the time for B° -> / (980)2Tg 
signal events. If an event has more than one candidate, 
we select one using a reproducible pseudo-random proce- 
dure based on the event timestamp. 

With the above selection criteria, we obtain a signal 
efficiency determined from Monte Carlo (MC) simulation 
of 21 — 25%, depending on the position in the DP. 

Of the selected signal events, 8% of B° -> p°K° s , 6% 
of B° -> if*(892) + 7r- and 4% of B° -> f (980)K% 
events are misreconstructed. Misreconstructed events 
occur when a track from the tagging B is assigned to 
the reconstructed signal candidate. This occurs most of- 
ten for low-momentum tracks and hence the misrecon- 
structed events are concentrated in the corners of the 
DP. Since these are also where the low-mass resonances 
overlap strongly with other resonances, it is important to 
model the misreconstructed events correctly. The model 
used to account for misreconstructed events is detailed 
in Sec. IVAl 

We use MC events to study the background from other 
B decays (B background). More than fifty channels were 
considered in preliminary studies, of which twenty are in- 
cluded in the final likelihood model - those with at least 
two events expected after selection. These exclusive B 



background modes are grouped into ten different classes 
that gather decays with similar kinematic and topological 
properties: nine for neutral B decays, one of which ac- 
counts for inclusive decays, and one for inclusive charged 
B decays. 

Table QT] summarizes the ten B background classes that 
are used in the fit. The yields of those classes that have 
a clear signature in the DP are allowed to float in the 
maximum-likelihood fit, the remainder are fixed. When 
the yield of a class is varied in the maximum-likelihood 
fit the quoted number of events corresponds to the fit 
results. For the other modes, the expected numbers of 
selected events are computed by multiplying the selection 
efficiencies (estimated using MC simulated decays) by the 
world average branching fractions [1, [3l| , scaled to the 
data set luminosity (347 fb _1 ). 



V. THE MAXIMUM-LIKELIHOOD FIT 



We perform an unbinned extended maximum- 
likelihood fit to extract the inclusive B° — > Kgir + ir~ 
event yield and the resonant amplitudes. The fit uses the 
variables toeSj AE, the NN output, and the SDP to dis- 
criminate signal from background. The At measurement 
allows the determination of mixing-induced CP violation 
and provides additional continuum background rejection. 

The selected on-resonance data sample is assumed to 
consist of signal, continuum background, and B back- 
ground components. The signal likelihood consists of the 
sum of a correctly reconstructed ("truth-matched," TM) 
term and a misreconstructed ("self-cross-feed," SCF) 
term. Generally, the components in the fit are separated 
by the flavor and tagging category of the tag side B de- 
cay. 

The probability density function (PDF) for an event 
i in tagging category c is the sum of the probability den- 
sities of all components, namely 



N- f c - 

J, sigJsig 



(1 - f SCF )K 



SCF) 1 sig-TM.i 



SCF ' sig-SCF.i 



+ N B+ f B+ -(l + qt a ^A B+ )V c B+ , t 



+ E *B° i fh j 'PZo# ■ (26) 

3=1 

The variables are defined in Table ED The PDFs 
V c x (X = {sig-TM, sig-SCF, qq, B+ , B }) are the 
product of the four PDFs of the discriminating variables, 
x\ = 7?t,es, x 2 = AE, X3 = NN output, and the triplet 
x 4 = {m',6',At}: 

4 

k=l 
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TABLE II: Summary of B background modes included in the fit model. When the yield is varied in the fit, the quoted number 
of events corresponds to the fit results. Otherwise, the expected number, taking into account the branching ratios and efficiency, 
is given. 







Vai*i oH 
V ell led 


BR 


l\Ti i m r»ov rif cnroYi^~c 
IN UlIlUcI Ul trvcllLb 


R° -> 

LJ r 


n~ ( > z<^o — \ 4- 


yes 




0077 _|_ flf) 


B° -» 


J/ty>(-> « + 0#s 


yes 




1803 ± 43 


B° - 




yes 




142 ± 13 


B° -» 




yes 




37 ± 16 


B° -> 




no 


(39.7 ±3.7) x 1(T 6 


7.3 ±0.7 


B° -» 


D7T)7T + 


no 


(2.57 ±0.10) x 10" 3 


43.8 ± 2.5 


B° -» 


D-h + ; 5° -> D-fi+Vp 


no 


(2.94 ±0.19) x 10~ 3 


281 ± 20 


B° — > 


D*-p+ 


no 


(14.2 ± 1.4) x 10" 3 


34.5 ±4.6 


B° — > 


{neutral generic decays} 


no 


not applicable 


114 ± 7 


B+ -» 


{charged generic decays} 


no 


not applicable 


282 ± 11 



TABLE III: Definitions of the different variables in the likelihood function given in Eq. (|26|l . 

Variable Definition 

iV s ig total number of K%ty + -k~ signal events in the data sample 

/si g fraction of signal events that are tagged in category c 

/scf fraction of SCF events in tagging category c, averaged over the DP 

^"sig-TM.i product of PDFs of the discriminating variables used in tagging category c for TM events 

"Psig-scF.i product of PDFs of the discriminating variables used in tagging category c for SCF events 

number of continuum events that are tagged in category c 

Qtag.i tag flavor of the event, defined to be ±1 for a -Bta g and —1 for a Bt ag 

A q q parameterizes possible asymmetry in continuum events 

Vqq ti continuum PDF for tagging category c 

-^ciass number of neutral _B-related background classes considered in the fit, namely nine 

N B+ number of expected charged B background events 

A^Oj number of expected events in the neutral B background class j 

fg + fraction of charged B background events that are tagged in category c 

fgOj fraction of neutral B background events of class j that are tagged in category c 

A B + describes a possible asymmetry in the charged B background 

V B+ i B + background PDF for tagging category c 

V B o y neutral B background PDF for tagging category c and class j 



where i is the event index and j is a B background class. 
Not all the PDFs depend on the tagging category; the 
general notations P£ ^ and V c x ^ are used for sim- 
plicity. Correlations between the tag and the position in 
the DP are absorbed in tag-flavor-dependent SDP PDFs 
that are used for continuum and charged B backgrounds. 
The parameters Ab+ and A q q parametrize any potential 
asymmetry between these PDFs. The extended likeli- 
hood over all tagging categories is given by 

7 N a 

£ = ii e ~* c n^> ( 2g ) 

c— 1 i 

where N c is the total number of events expected in cat- 
egory c. 

A total of 75 parameters are varied in the fit. They in- 
clude the 12 inclusive yields (signal, four B background 



classes, and seven continuum yields, one per tagging cat- 
egory), 30 parameters for the complex amplitudes from 
Eq. and 33 parameters of the different PDFs. The 

latter include most of the parameters describing the con- 
tinuum distributions. 



A. The At and Dalitz plot PDFs 

The SDP PDFs require as input the DP-dependent 
selection efficiency, e — e(m',8'), and SCF fraction, 
/scf = /scf(w', #'). Both quantities are taken from 
MC simulation. Away from the DP corners the efficiency 
is uniform. It decreases when approaching the corners, 
where one of the three particles in the final state is nearly 
at rest so that the acceptance requirements on the par- 
ticle reconstruction become restrictive. Combinatorial 
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backgrounds and hence SCF fractions are large in the 
corners of the DP due to the presence of soft tracks. 
For an event i we define the time-dependent SDP PDFs 

P sig _TM,i(m',0',Ai)= (29) 
e l (l-/scF, l )|detJ l ||^ ± (A0| 2 , 

P B ig-scF,i(m , ,0 / ,At)= (30) 
^/scF^ldetJ.II^At)! 2 , 

where P s ig-TM,i{m' , 9', At) and P sig -scF,i(m' , 9', At) are 
normalized to unity. The phase space integration involves 
the expectation values (e(l — /scf) | det J\ FkF£,) and 
(e /scf I det J\ F k Fj*,) for TM and SCF events, where the 
indices k, k! run over all resonances belonging to the sig- 
nal model. The expectation values are model-dependent 
and are computed by MC integration over the SDP: 

(£(l-/scF)|detJ|F feJ F fe *,) = (31) 
Jo Jo £ U ~ /scf) I det J\ F k F£, dm'd9' 
/* \l e | det J| F fe F* dmW 

and similarly for (e/scF | det J\ F k F£,), where all quan- 
tities in the integrands are DP-dependcnt. 

Equation (f2"fi)) invokes the phase space-averaged SCF 
fraction / SCF = (/ SCF | det J\ F k F*,). The PDF normal- 
ization is decay-dynamics-dependent and is computed it- 
eratively. We determine the average SCF fractions sepa- 
rately for each tagging category from MC simulation. 

The width of the dominant resonances are large com- 
pared to the mass resolution for TM events (about 
8MeV/c 2 core Gaussian resolution). We therefore neglect 
resolution effects in the TM model. Misreconstructed 
events have a poor mass resolution that strongly varies 
across the DP. It is described in the fit by a 2 x 2- 
dimensional resolution function 

R SCF (m' rl 8' r ,m' t ,0' t ) , (32) 

which represents the probability to reconstruct at the co- 
ordinates (m' r , 6' r ) an event that has the true coordinates 
(m' t ,9' t ). It obeys the unitarity condition 

l l 

R S CF(m' r ,6' r ,m' u 9' t )dm' r d0' r = 1 , (33) 
o o 

and is convolved with the signal model. The -Rscf func- 
tion is obtained from MC simulation. 

We use the signal model described in Sec. Ill Al It con- 
tains the dynamical information and is connected with 
At via the matrix element in Eq. (|22p . which intervenes 
in the signal PDFs defined in Eq. ^ and §U§. The 
PDFs arc diluted by the effects of mistagging and the 
limited vertex resolution (42j. The At resolution func- 
tion for signal (both TM and SCF) and B background 
events is a sum of three Gaussian distributions. The pa- 
rameters of the signal resolution function are determined 
by a fit to fully reconstructed B° decays [40j |. 



The charged B background contribution to the likeli- 
hood, given in Eq. (f2"rJ|) , uses distinct SDP PDFs for each 
reconstructed B flavor tag, and a flavor-tag-averaged 
PDF for untagged events. The PDFs are obtained from 
MC simulation and are described by histograms. The 
At resolution parameters are determined by a fit to fully 
reconstructed B + decays. For the B + background class 
we adjust the effective lifetime to account for the misre- 
construction of the event that modifies the nominal At 
resolution function. 

The neutral B background is parameterized with PDFs 
that depend on the flavor tag of the event. In the case of 
CP eigenstates, correlations between the flavor tag and 
the Dalitz coordinates are expected to be small. How- 
ever, non-CP eigenstates, such as afir T , may exhibit 
such correlations. Both types of decays can have di- 
rect and mixing-induced CP violation. A third type of 
decay involves charged D mesons and does not exhibit 
mixing-induced CP violation, but usually has a strong 
correlation between the flavor tag and the DP coordi- 
nates because it consists of £?-flavor eigenstates. Direct 
CP violation is also possible in these decays, though it 
is set to zero in the nominal model. The DP PDFs are 
obtained from MC simulation and are described by his- 
tograms. For neutral B background, the signal At res- 
olution model is assumed. Note that the SDP- and At- 
dependcnt PDFs factorize for the charged B background 
modes, but not necessarily for the neutral B background 
due to B°B° mixing. 

The DP treatment of the continuum events is similar 
to that used for charged B background. The SDP PDF 
for continuum background is obtained from on-resonance 
events selected in the Toes sidebands and corrected for 
feed-through from B decays. A large number of cross 
checks have been performed to validate the empirical 
shape used. The continuum At distribution is param- 
eterized as the sum of three Gaussian distributions with 
common mean and three distinct widths that scale with 
the At per-event error. This introduces six shape pa- 
rameters that are determined by the fit. The model is 
motivated by the observation that the mean of the At 
distribution is independent of the per-event error, and 
that the width depends linearly on this error. 



B. Description of the other variables 

The toes distribution of TM signal events is parame- 
terized by a bifurcated Crystal Ball function [43I.I44L l45j ] , 
which is a combination of a one-sided Gaussian and a 
Crystal Ball function. The mean and the two widths of 
this function are determined by the fit. The AE distri- 
bution of TM signal events is parameterized by a double 
Gaussian function. The five parameters of this function 
are determined by the fit. Both toes an d AE PDFs 
are described by histograms, taken from the distributions 
found in appropriate MC samples, for SCF signal events 
and all B background classes. Exceptions to this are the 
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FIG. 1: Standard (left) and square (right) Dalitz plots of the selected data sample of 22525 B — > K%-k + -k candidates. The 
narrow bands correspond to D^ 1 ^, J/'ipKg and ip(2S)Kg background events. 



m E s PDFs for the B° -> D^7r+ and B° -> J/0 -Kg com- 
ponents, and the AE PDF for B° — ► Z?~7r + , which are 
the same as the corresponding distributions of TM signal 
events. The toes and AE PDFs for continuum events are 
parameterized by an ARGUS shape function [46[ and a 
first-order polynomial, respectively, with parameters de- 
termined by the fit. 

We use histograms to empirically describe the distri- 
butions of the NN output found in the MC simulation 
for TM and SCF signal events and for all B background 
classes. We distinguish tagging categories for TM signal 
events to account for differences observed in the shapes. 
The continuum NN distribution is parameterized by a 
third-order polynomial that is constrained to take pos- 
itive values in the range populated by the data. The 
coefficients of the polynomial are determined by the fit. 
Continuum events exhibit a correlation between the DP 
coordinates and the shape of the event that is exploited 
in the NN. To correct for this effect, we introduce a 
linear dependence of the polynomial coefficients on the 
variable Adp, defined as the smallest of the three in- 
variant masses, and is thus a measure of the distance 
of the DP coordinates from the kinematic boundaries of 
the DP. The parameters describing this dependence are 
determined by the fit. 



VI. FIT RESULTS 

The standard and square Dalitz plots of the se- 
lected data sample are shown in Fig. [1] The 
maximum-likelihood fit of 22525 candidates results in a 
B° — > Kgir + ir~ event yield of 2182±64 and a continuum 
yield of 14240 ± 126, where the uncertainties are statisti- 
cal only. The remaining number of events is covered by 
the yields of backgrounds from charged and neutral B 
decays, where the dominant contributions are 3361 ± 60 



B° -> L>"tt+ and 1804 ± 44 B° -> J/ipK% events. 

When the fit is repeated starting from input parame- 
ter values randomly chosen within wide ranges above and 
below the nominal values for the magnitudes and within 
the [0 — 360°] interval for the phases, we observe conver- 
gence toward two solutions with minimum values of the 
negative log likelihood function — 21og£ that are equal 
within 0.32 units. In the following, we refer to them as 
solution I (the global minimum) and solution II (a local 
minimum). Between the two solutions, the fit values for 
most free parameters are very similar. Exceptions occur 
among isobar parameters, and most particularly isobar 
phases, some of which can differ significantly. 

For a given event i, we define the likelihood ratio as 
R = T sig-TM, i/'Pi (see Eq. (|2"o) and explanations below). 
Figure [2] shows distributions of log R for all the events en- 
tering the fit, and for the signal-like region. We obtain 
signal enriched samples that are used in some of the fig- 
ures below, by removing events with small values of R; 
in each case R is computed excluding the variable being 
plotted. Figure [3] shows distributions of AE, toes, an d 
the NN output which are enhanced in signal content by 
requirements on R. Figures [4] to [7] show similar distribu- 
tions for TO(7r + 7r _ ), m(K®ir), and Adp- These distribu- 
tions illustrate the good quality of the fit in the signal- 
enhanced regions. Signal enriched distributions of At and 
At asymmetry for events in the regions of /o (980) Kg and 
p°(770)K^ are shown in Fig. [SJ 

In the fit, we measure directly the relative magni- 
tudes and phases of the different components of the sig- 
nal model. The magnitude and phase of the B° — > 
/o(980).Kly amplitude are fixed to 4 and 0, respectively, 
as a reference. The results corresponding to the two so- 
lutions are given together with their statistical uncer- 
tainties in Table [TVl The full (statistical, systematic and 
model dependent) correlation matrices between the mag- 
nitudes and the phases for the two solutions are given 
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FIG. 2: Distributions of the logarithm of likelihood ratio (log R) for all events entering the fit (left) and in the signal-like 
region (right). In the right hand side plot, a veto in the D~rt + , J/tpKg, and ip(2S)K < g bands has been applied. Points with 
error bars give the on-resonance data. The solid histogram shows the projection of the fit result. The dark, medium, and 
light shaded areas represent respectively the contribution from continuum events, the sum of continuum events and the B 
background expectation, and the sum of these and the misreconstructed signal events. The last contribution is hardly visible 
due to its small fraction. Below each bin are shown the residuals, normalized in error units. The parallel dotted and full lines 
are the 1<t and 2<r deviations. Points, histograms, shaded areas, and residual plots have similar definitions in Fig. |3]to|8] 



TABLE IV: Results of fit to data for the isobar amplitudes with statistical uncertainties. Both solutions are shown. 









Solut 


ion I 






Solution II 


Isobar Amplitude Magnitude 


Phase (°) 


Mag 


nitude 


Phase (°) 


C /o(980)if° 




4.0 


0.0 




4.0 


0.0 


C /o(980)Kg 


3.7 


± 


0.4 


-73.9 ± 19.6 


3.2 


± 


0.6 


-112.3 ± 20.9 


C p(770)K° 


0.10 


± 


0.02 


35.6 ± 14.9 


0.09 


± 


0.02 


66.7 ± 18.3 


Cp(770)Kg 


0.11 


± 


0.02 


15.3 ±20.0 


0.10 


± 


0.03 


-0.1 ± 18.2 


c X*+(892)ir- 


0.154 


± 


0.016 


-138.7 ±25.7 


0.145 


± 


0.017 


-107.0 ±24.1 


C K *- (892)vr+ 


0.125 


± 


0.015 


163.1 ±23.0 


0.119 


± 


0.015 


76.4 ± 23.0 




6.9 


± 


0.6 


-151.7 ± 19.7 


6.5 


± 


0.6 


-122.5 ± 20.3 




7.6 


± 


0.6 


136.2 ± 19.8 


7.3 


± 


0.7 


52.6 ±20.3 


c / 2 (1270)KTg 


0.014 


± 


0.002 


5.8 ± 19.2 


0.012 


± 


0.003 


23.9 ±22.7 


C/ 2 (1270)if« 


0.011 


± 


0.003 


-24.0 ±28.0 


0.011 


± 


0.003 


-83.3 ± 24.3 


c /x(1300)Kg 


1.41 


± 


0.23 


43.2 ± 22.0 


1.40 


± 


0.28 


85.9 ± 24.8 


C /x(1300)JfO 


1.24 


± 


0.27 


31.6 ±23.0 


1.02 


± 


0.33 


-67.9 ± 22.1 


CNR 


2.6 


± 


0.5 


35.3 ± 16.4 


1.9 


± 


0.7 


56.7 ± 23.6 


CNR 


2.7 


± 


0.6 


36.1 ±18.3 


3.1 


± 


0.6 


-45.2 ± 17.8 




0.33 


± 


0.15 


61.4 ± 44.5 


0.28 


± 


0.16 


51.9 ±38.4 


c x. a i<% 


0.44 


± 


0.09 


15.1 ±30.0 


0.43 


± 


0.08 


-58.5 ±27.9 



in the Appendix. The measured relative amplitudes Cfe, fined as: 

where the index represents an intermediate resonance, 2 

are used to extract the Q2B parameters denned below. C(k) = ^ Cfc ^ ~ ^ Cfc ^ (35) 

For a resonant decay mode k which is a CP eigenstate, l c fe| 2 + \^k\ 

the following Q2B parameters are extracted: the angle 
/3 e ff defined as 



/3 eff (fc) = iarg( C)t c£) (34) 



S(k) = . (36) 



For a flavor-specific resonant decay mode k such as 
and the direct and mixing- induced CP asymmetries, de- B° — > A"* + (892)7r~, it is customary to define the direct 
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FIG. 3: Distributions of AE (left), w,es (center), and AW output (right) for a sample enhanced in B° — > A"57r + 7r~ signal with 
a requirement on the likelihood ratio R computed without the variable being plotted. In each case the applied cut rejects 99% 
of continuum background, while retaining 28% of signal for AE and jtibs, and 16% for NN. A veto in the D~tt + and J/%j)K% 
bands has been applied. 




FIG. 4: Spectra of m^+^.- (left) and symmetrized m K o w (right) for the whole data sample. For m^+^-, the insets show the 
J/tp region (a) and in the ip(2S) region (b). The symmetrized rn K ^ is obtained by folding the SDP with respect to the 6' 
variable at 0.5. The inset shows the D region. 




FIG. 5: Distribution of m 7r + 7r - for a sample enhanced in B° — > Kg-K + TT signal, showing the /o (980) .Kg and p°{77Q)K% signal 
region for positive (left) and negative (right) 7r + 7r~ helicity. The contribution from f x {1300)K° s and / 2 (1270)A"g are also 
visible. A veto in the Dn band has been applied. The At and DP PDFs have been excluded from the likelihood ratio R used 
to enhance the sample in signal events. The cut on R retains 21% of signal, while rejecting 99% of continuum. 
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FIG. 6: Distributions of rn^^ for a sample enhanced in B° — ► Kgir + Tv signal, showing the K*(892)iv and _K"*(1430)7r signal 

region for positive (left) and negative (right) K%n helicity. A veto in the J/t/jKg and ip{2S)K% bands has been applied. The 
At and DP PDFs have been excluded from the definition of the likelihood ratio used to enhance the sample in signal events. 
The cut on R retains 18% of signal while rejecting 94% of continuum. An interference between the vector and scalar K* + is 
apparent through a positive (negative) forward-backward asymmetry below (above) the .ft"* (892). 
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FIG. 7: Distributions of the Adp variable, for a sample enhanced in B° — > Kgir + Tr~ signal. The Adp variable is defined as 
min(m K o n + , rn K o n - , ra 7T + 7V - ). Small (large) values of Adp correspond to the edges (center) of the DP. On the left (right) side 

of the figure, for Adp < 1.9GeV/c 2 (> 1.9GeV/c 2 ), the dominant contribution to the signal is from the light resonances (the 
NR) component of the signal model. A veto in the Dir, J/ipKg, and ip(2S)K% bands has been applied. The At and DP PDFs 
have been excluded from the likelihood ratio R used to enhance the sample in signal events. The cut on R retains 37% of signal 
while rejecting 88% of continuum. 
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FIG. 8: Distributions of At when the B^ is a B° (top), B° (middle), and the derived At asymmetry (bottom). Plots on the 
left (right) hand side, correspond to events in the fo(980)Kg (p°(770)K s ) region. These distributions correspond to samples 
where the D~ir + and J/ip Kg bands are removed from the DP, and the At and DP PDFs have been excluded from the likelihood 
ratio R used to enhance the sample in signal events. The cut on R retains 24% of signal while rejecting 98% of continuum. 



CP asymmetry parameter Acp as 
A C p(k) 



c- k \ 2 ~\c k \ 2 



+ \c k \ 



(37) 



For a pair of resonances k and k', the phase 4>(k, k') 
relating their amplitudes Ck and Ck> , defined as 



4>{k, k') = arg(cfeCfc,), 



(38) 



can be accessed by exploiting the interference pattern in 
the DP areas where k and k! overlap; correspondingly, 
the phase 4>(k, k') for the CP-conjugated amplitudes Ck 
and Ck' is 



<t>{k,k') = arg(c fc 4,). 



From these two phases, the difference A0(fc, k') = 
4>{k, k') — 4>(k, k'), can be extracted. This parameter is a 



direct CP violation observable, and can only be accessed 
in an amplitude analysis. 

For a resonant decay mode k, the phase relating its 
amplitude Ck to its charge conjugate Ck is defined as 



A$(fc) = arg(c fc c£); 



(40) 



here it is worth recalling that we use a convention in 
which the B° decay amplitudes have absorbed the phase 
from B°B° mixing, and so the phase of q/p is implicit in 
the A§(k) parameter. Although the definition of this 
parameter is technically similar to the /3 e ff phase de- 
fined in Eq. (|34[) . they differ in their physical interpreta- 
tion. The parameter /3 e ff quantifies the time-dependent 
mixing-induced CP asymmetry, and therefore is most rel- 



(39) evant for the CP eigenstate modes, such as p°(770)Kg 



and / (980)ifg. On the other hand the A$(fc) pa- 
rameter concerns mostly flavor-specific modes, such as 
B° — > i^* + (892)7r _ , for which there is no interference be- 
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tween decays with and without mixing. For such modes, 
sensitivity to A$(fc) is provided indirectly by the inter- 
ference pattern of the resonance k with other modes that 
are accessible both to B° and B° decays. 

We also extract the relative fit fraction FF of a Q2B 
channel k, which is calculated as: 



FFik) = 



c k 



Euv (c^ + ^c*)^*) 



where the terms 



F F 



'ds+ds- 



(41) 



(42) 



are obtained by integration over the complete Dalitz plot. 
The total fit fraction is defined as the algebraic sum of all 
fit fractions. This quantity is not necessarily unity due to 
the potential presence of net constructive or destructive 
interference. Using the relative fit fractions, we calculate 
the branching fraction B for the intermediate mode k as 



FF(k) x B(B° -> K a s n+n-) 



(43) 



where B(B° — > K s tt + tt ) is the total inclusive branching 
fraction 



B(B° 



(44) 



BB 



We compute the average efficiency, e, by weighting MC 
events with the measured intensity distribution of sig- 
nal events, (|.4(DP)| 2 + |„4_(DP)| 2 )/2. The term N BB 
is the total number of BB pairs in the sample. Fi- 
nally, we use the following integrals of amplitudes over 
the complete Dalitz plot to measure the inclusive direct 
CP-asymmetry: 



a incl 
A CP 



Jf{\ 


A\ 


2 - I-4I 2 ) 


1 ds+ds- 


//(I 


A\ 


2 + l< 


| ds + ds_ 



(45) 



The Q2B parameters and fit fractions are given in Ta- 
ble [V] together with their statistical and systematic er- 
rors. The branching fractions are shown in Table IVT1 

To extract the statistical uncertainties on the Q2B pa- 
rameters we perform likelihood scans, not relying on any 
assumption about the shape of the likelihood function. 
Since the Q2B parameters are not directly used in the fit, 
we instead must perform the scan fixing one or two pa- 
rameters among the signal model magnitudes and phases. 
These are chosen in such a way that the resulting like- 
lihood curve can be trivially interpreted in terms of the 
Q2B parameter of interest. In each case the chosen pa- 
rameters are fixed at several consecutive values, for each 
of which the fit to the data is repeated. The error on 
the Q2B parameter is determined by the points, or the 
contour, where the —2 log£ function changes by one unit 
with respect to its minimum value. Systematic uncertain- 
ties are discussed in Sec. IVII1 Results of the likelihood 
scans in terms of — 2 A log £ are shown in Fig. [91 to [TBI 



The measurements of time-dependent CP-violation in 
the f (980)K s and p (770)K° s modes are presented as 
two-dimensional likelihood scans in the (/3 e ff> C) plane, 
shown in Fig. [5] The scans are displayed as confidence 
level contours after two-dimensional convolution with the 
covariance matrix of systematic uncertainties. On the 
same figure are also displayed the one-dimensional like- 
lihood scans of (3 e g. For f (980)K s the two solutions lie 
below and above 45 degrees and correspond very closely 
to the trigonometric ambiguity between a given value of 
(3 e g and 90° — (3 c g (mirror solutions). On the other hand, 
for p° (770) Kg both solutions are below 45 degrees. In 
this case the local solutions corresponding to the trigono- 
metric ambiguities of the two observed solutions are sup- 
pressed at 3.6 and 2.0 standard deviations, respectively. 

The ((3 e ff,C) plane can be transformed to the more 
familiar (5, C) plane using Eq. (|3"4")) to (f3"6")l . The corre- 
sponding two-dimensional contours are shown in Fig. 1101 
While a part of the information on the phases is lost, this 
representation has nonetheless the advantage of allowing 
direct comparison with the measurement of sin 2(3 and C 
in b — ► ccs modes. For fo(980)K s , the results agree with 
the expectation based on b — ► ccs to I.Ict; for p° (770) Kg 
the agreement is better than la. For the measured values 
of ((3 e H,C) for f (980)K s , CP conservation is excluded 
at 3.5ct. For p°(770)K s , the measurement of ((3 C «, C) is 
consistent with CP conservation within la. 

The measurement of the phase A$(K* + (892)n~) is 
presented as a one-dimensional likelihood scan in Fig. 1111 
For this flavor-specific mode, there is virtually no region 
in phase space that is accessible both to B° and B°; thus, 
sensitivity to this phase difference is limited. Simula- 
tion shows that interference of the K* + (892)ir~ with the 
fo(980)K s and p°(770)K s modes (for which B° and B° 
amplitudes interfere via mixing) provides most of the sen- 
sitivity to A$(if*+ (892)7r~); unfortunately, the overlap 
in phase space of these resonances is small. As a conse- 
quence, only the (—137, —5)° interval is excluded at 95% 
confidence level. Figure [Til also shows the measurement 
of the similar phase difference for the (Ktt)q component. 
As for K*(892), the measurement sets no strong con- 
straint on this phase. Only the interval [—132, +25]° is 
excluded at 95% confidence level. 

In contrast, due to the sizable overlap in phase 
space between the Kir S- and P- waves of the same 
charge, the relative phases </)((Ktt)q , K* ± (892)) are 
measured to ±13° including systematics. The one- 
dimensional scans are shown in Fig. [12J The associated 
observable A^Kir)^ , K* ± (892)) is compatible with 
CP conservation. Figure [T2] also shows the scans for 
#/o(980)tfg, p°(770)K° s ), c/>(p (770)K° s , K*+ (892)^) , 
and their corresponding CP-conjugates. It is clear from 
this figure and from Table V that the phases for the for- 
mer are measured to a better accuracy. This is due to 
the larger overlap in phase space between the /o(980) and 
the p°(770). In both cases the associated observables A(j) 
are compatible with CP conservation. 

For the remaining resonant modes in the signal DP 
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TABLE V: Summary of measurements of the Q2B parameters for solutions I and II. The first uncertainty is statistical, 
the second is systematic, and the third represents the DP signal model dependence. We also show the total (statistical and 
systematic) linear correlations between the parameters f3 c g (S) and C. Phases are given in degrees and FFs in percent. 



Parameter 


Solution I 


Solution II 


C(/ (980)X°) 
/3eff(/o(980)^§) 
S(f (980) Kg) 

CoTT[/3 eS (M980)K s ),C(f (980)K s )] 

Corr[S(/o(980)#g),C(/ (980)Xg)] 

FF(f (980)K° s ) 


0.08 ±0.19 ±0.03 ±0.04 

36.0 ±9.8 ± 2.1 ± 2.1 
-0.96+Q 04 ± 0.03 ± 0.02 
-3.1% 
19.7% 
13.8±i t ±0.8 ±0.6 


0.23 ±0.19 ±0.03 ±0.04 
56. 2 ±10. 4 ±2.1 ±2.1 
-0.90lo o8 ± 0.03 ± 0.02 
'-17.0% 
12.5% 
13.5±i'3 ±0.8 ±0.6 


C(p°(770)K° s ) 

/3eff(p°(770)7fg) 

S(p°(770)K° s ) 

Corr[/3 eS (p (770)K s ),C(p (770)K° s )] 
Cott[S(p (77Q)K s ),C(p (77Q)K° s )] 
FF(p° (770) K° s ) 


-0.05 ±0.26 ±0.10 ±0.03 
10.2 ±8.9 ±3.0 ±1.9 
0.35tg;^ ± 0.06 ± 0.03 
-23.0% 
-21.3% 
8.6li;3 ±0.5 ±0.2 


-0.14 ±0.26 ±0.10 ±0.03 
33.4 ±10.4 ±3.0 ±1.9 
0.9ltg^ ± 0.06 ± 0.03 
-34.0% 
-10.4% 
8.5±i;l ± 0.5 ± 0.2 


A CP (K* (892)n) 

A$(7f*(892)7r) 

FF(K*(892)tt) 


-0.21 ±0.10 ±0.01 ±0.02 
58.3 ±32.7 ±4.6 ±8.1 
n.otl l ±0.6 ±0.8 


-0.19^ 11 ±0.01 ±0.02 
176.6 ± 28.8 ±4.6 ±8.1 
10.9ti'o ±0.6 ±0.8 


A CP {(Ktt)otv) 
FF((Kv)fa) 


0.09 ±0.07 ±0.02 ±0.02 
72.2 ±24.6 ±4.1 ±4.4 
45. 2 ±2. 3 ±1.9 ±0.9 


0.12+°-^ ±0.02 ±0.02 
-175.1 ±22.6 ±4.1 ±4.4 
46.1 ±2.4 ±1.9 ±0.9 



C(/ 2 (1270)Kg) 
/3 eff (/ 2 (1270)7f§) 


0.28±g-^ ± 0.08 ± 0.07 


0.09 ±0.46 ±0.08 ±0.07 


14.9 ±17.9 ±3.1 ±5.2 


53.6 ±16.7 ±3.1 ±5.2 


S(M1270)K° S ) 


-0.48 ±0.52 ±0.06 ±0.10 


-0.95 ±0.17 ±0.06 ±0.10 


Corr[/3 eff (/ 2 (1270)^),C(/ 2 (1270)^§)] 


11.5% 


-2.8% 


Corr[5(/ 2 (1270)^§),C(/ 2 (1270)^§)] 


0.9% 


21.2% 


FF(f 2 (l270)K° s ) 


2.3t 1 ;? ± 0.2 ± 0.7 


2.3t[]; l 7 ± 0.2 ± 0.7 



C(fx(1300)K° s ) 


0.13t°'35 ±0.04 ±0.09 


0.30t^t ± 0- 04 ± O- 09 


l3 cH (fx(1300)K° s ) 


5.8 ±15.2 ±2.2 ±2.3 


76.9 ±13.8 ±2.2 ±2.3 


S(fx(1300)K° s ) 


-0.20 ±0.52 ±0.07 ±0.07 


-0.42 ±0.41 ±0.07 ±0.07 


Corr [f5 eB (f x ( 1300) tfg ),C(fx (1300) Jfg)] 


-27.0% 


-9.3% 


Corr[S(/x(1300)Kg),C(/x(1300)^g)] 


28.5% 


6.1% 


FF(f x (1300)Xg) 


3.6t l ± 0.3 ± 0.9 


3.5+Vg ± 0.3 ± 0.9 



C(NR) 
Pes (NR) 
S{NR) 

CoTT[/3 eS {NR),C(NR)] 

Corr[S(NR),C(NR)] 

FF(NR) 


0.01 ±0.25 ±0.06 ±0.05 

0.4 ± 8.8 ± 1.9 ±3.8 
-0.01 ±0.31 ±0.05 ±0.09 
-10.6% 
10.6% 
11. 5 ±2.0 ±1.0 ±0.6 


-0.45^24 ± 0.06 ± 0.05 
51.0 ±13.3 ±1.9 ±3.8 
-0.87 ±0.18 ±0.05 ±0.09 
-37.9% 
-91.5% 
12. 6 ±2.0 ±1.0 ±0.6 


C( X coK° s ) 

Pes(XcoK° s ) 

S(XcoK° s ) 

CoTT[/3 eS (xc K s ),C(xc K s )} 

Cott[S(xcoK° s ),C(xc K s )] 

FF(xc K° s ) 


-0.29±o44 ± 0.03 ± 0.05 
23.2 ±22.4 ±2.3 ±4.2 
-0.69 ±0.52 ±0.04 ±0.07 
-5.8% 
-19.1% 
IMto fs ±0.04 ±0.11 


-0.41±o1 2 ±0.03 ±0.05 
55.2 ±23.3 ±2.3 ±4.2 
-0.85 ±0.34 ±0.04 ±0.07 
-5.8% 
-74.2% 
0.99l°;3o ±0.04 ±0.11 


total FF 

a incl 
A CP 


97.2^-^2.1 ± 1.15 
-0.01 ±0^05 ±0.01 ±0.01 


98.3±1;| ±2.1 ± 1.15 
0.01 ±0.05 ±0.01 ±0.01 


4>(J° (980)K° S , p(77Q)K° s ) 
<t>{K* {892)tt, (Kny 7v) 
ct>{p{770)K%,K*{892)^) 
4>(p(770)K° s , (Kny o n) 


-35.6 ±14.9 ±6.1 ±4.4 
13.0 ± 10.9 ±4.6 ±4.7 
174.3 ±28.0 ±8.7 ±12.7 
-172.8 ±22.6 ±10.1 ±8.7 


-66.7 ±18.3 ±6.1 ±4.4 
15.5 ± 10.2 ±4.6 ±4.7 
-173.7 ± 29.8 ± 8.7 ± 12.7 
-170.8 ±26.8 ±10.1 ±8.7 


t(f(980)K s ,p(770)K° s ) 
4>(K* (892)tt, (Ktt)^) 
0(p(77O)K° s ,K*(892)iv) 
0(p(77O)X§, (A-jtJStt) 


-89.2 ±17.1 ±8.5 ±7.2 
26.9 ±9.2 ±4.9 ±6.1 
-147.8 ± 24.7 ± 11.3 ± 11.9 
-120.9 ±21.6 ±8.7 ±7.3 


-112.2 ±17.8 ±8.5 ±7.2 

23.8 ±9.1 ±4.9 ±6.1 
-76.5 ± 24.0 ± 11.3 ± 11.9 
-52.7 ±21.4 ±8.7 ±7.3 
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TABLE VI: Summary of measurements of branching fractions averaged over charge conjugate states. The quoted numbers were 
obtained by multiplying the corresponding fit fractions by the measured inclusive B° — > K°tt~it~ branching fraction. R denotes 
an intermediate resonant state and h stands for a final state hadron: a charged pion or a K°. To correct for the secondary 
branching fractions we used the values from Ref. [31| and B(K* + (892) — > K°ir + ) — |. The first uncertainty is statistical, the 
second is systematic, and the third represents the DP signal model dependence. The fourth errors, when applicable, are due 
to the uncertainties on the secondary branching fractions. The quoted central values correspond to the global minimum, and 
errors account for the presence of the second solution. 



Mode 


B(B° -» Mode) x B(R -> hh) x 10" 6 


B(B° Mode) x 10" 6 


Inclusive B a — > K°ir + ir~ 




50.15 ±1.47 ±1.60 ±0.73 


fo{980)K° 


6.92 ± 0.77 ± 0.46 ± 0.32 




p"(770)K° 


4.311^6? ± 0.29 ± 0.12 


4.36±{J;£a ± 0.29 ± 0.12 ± 0.01 


K*+(892)ty- 


5.52lo;s4 ± 0.35 ± 0.41 


8.29l°;gi ± 0.53 ± 0.62 




22.7li;3 ± 1-2 ±0.6 




/ 2 (1270) if 


1.15±g;|| ±0.11 ±0.35 


2.7ilg:S±o.26±o.83tg:8S 


fx(1300)K° 


1.8llg;^ ±0.16 ±0.45 




flat NR 




5 - 77 ii:oo ±0.53 ±0.31 


XcoK 


0.52lo:?g ± 0.03 ± 0.06 


U2+H ± 8 ± 16 ± 12 
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model: f x (1300)K%, f 2 {l27Q)K% NR, and X co(lP)K° s , 
we scan the likelihood as a function of the corresponding 
fit fractions. These scans are shown in Fig. [131 We ob- 
tain a total (statistical and systematic) significance of 4.8 
and 3.8 standard deviations for the NR and XcoO-P)Kg 
components, respectively. The significance for the sum 
of fit fractions of the / 2 (1270) J FsTg and / x (1300)Kg com- 
ponents is 4.8 standard deviations while their individual 
significances are 2.9a and 2.4c, respectively. 

The (Ktt)q component is modeled in our analysis by 
the LASS parametrization (34J, which consists of a NR 
effective range term plus a relativistic Breit-Wigner term 
for the if* (1430) resonance. We separate from the 
corresponding branching fraction, quoted in Table IVI1 
the contribution of the K*(1A30) resonance and find it 
to be (29.9±^ ± 1.6 ± 0.6 ± 3.2) x 10" 6 . This value 
is corrected for the secondary branching fraction using 
B{K*(U30) -> Kir) from Ref. [H[ and the isospin rela- 
tion B(K*+(U30) -> K°tt+)/B(K*+(U30) -> K+n ) = 
2. The first uncertainty is statistical, the second is sys- 
tematic, the third represents the DP signal model de- 
pendence, and the fourth is due to the uncertainty on 
the secondary branching fraction. In addition we calcu- 
late the total NR contribution by combining coherently 
the effective range part of the LASS parametrization and 
the flat phase-space NR component. We find this total 
NR fit fraction to be 22.lt|§ ± 2.1 ± 0.7%. Note that 
this number accounts for the destructive interference be- 
tween the two NR terms. The corresponding branching 
fraction is (11.07±ojg ± 0.81 ± 0.40) x 10~ 6 . 

As a validation of our treatment of the time- 
dependence, we allow Tgo and Amj to vary in the fit. We 
findr B o = 1.579±0.061psand Am d = 0.497±0.035 ps" 1 
while the remaining free parameters are consistent with 
the nominal fit. The numbers for Tgo and Amd are 
in agreement with current world averages Q. In ad- 
dition we perform a fit floating the S parameters for 
5° -> J/i>K% and B° -> $(2S)K% events. We find 
S = sin(2/3) = 0.690 ± 0.077 and 0.73 ± 0.27 for J^K% 
and ip(2S)Kg respectively. These numbers are in agree- 
ment with the current world average [|| . Signal enhanced 
distributions of At and the At asymmetry for events in 
the J/ipKg region are shown in Fig. [T7] To validate the 
SCF modeling, we leave the average SCF fractions per 
tagging category free to vary in the fit and find results 
that are consistent with the MC estimation. 

As a further cross-check of the results, we performed 
an independent analysis and obtained compatible re- 
sults [47J • The main differences between this cross-check 
analysis and the one presented here were the use of 
a Fisher discriminant instead of a NN, the removal of 
bands in invariant mass to cut away the B — > D~ir + , 
B° -> J/ipK% and B° -> 4>(2S)K% contributions, and 
the use of Cartesian isobar parameters. 



VII. SYSTEMATIC STUDIES 

To estimate the contribution to B° — > K^tt + ii~ decay 
via other resonances, we first fit the data including these 
other decays in the fit model. We consider possible res- 
onances, including w(782), p°(1450), p°(1700), /o(1710), 
/ 2 (1810), ^±(1680), 2^(1430), X&QP), and a low 
mass a. A relativistic Breit-Wigner lineshape is used to 
parameterize these additional resonances, with masses 
and widths from Ref. [3lJ. As a second step we simu- 
late high statistic samples of events, using a model based 
on the previous fits, including the additional resonances. 
Finally, we fit these simulated samples using the nom- 
inal signal model. The systematic effect (contained in 
the "DP Model" field in Table |VIlJ) is taken from the 
difference observed between the generated and fitted val- 
ues. We quote this DP model uncertainty separately from 
other systematics. 

We vary the mass, width, and any other parameters of 
all isobar fit components within their errors, as quoted 
m Table HI and assign the observed differences in the 
measured amplitudes as systematic uncertainties ( "Line- 
shape" in Table ELD). 

To validate the fitting tool, we perform fits on large 
MC samples of fully-reconstructed events with the mea- 
sured proportions of signal, continuum, and B back- 
ground events. No significant biases are observed in these 
fits and therefore no corrections are applied. The sta- 
tistical uncertainties on the fit parameters are taken as 
systematic uncertainties ("Fit Bias" in Table IVTI1) . 

Another major source of systematic uncertainty is the 
B background model. The expected event yields from 
the background modes are varied according to the uncer- 
tainties in the measured or estimated branching fractions. 
Since B background modes may exhibit CP violation, the 
corresponding parameters are varied within their uncer- 
tainties, or, if unknown, within the physical range. As is 
done for the signal PDFs, we vary the At resolution pa- 
rameters and the flavor-tagging parameters within their 
uncertainties and assign the differences observed in these 
fits with respect to the nominal fit as systematic errors. 
These errors are listed as U B Background " in Table lyTIl 

Other systematic effects are much less important for 
the measurements of the amplitudes and are combined in 
the "Other" field in Table IVIII Details are given below. 

The parameters of the continuum PDFs are deter- 
mined by the fit. No additional systematic uncertain- 
ties are assigned to them. An exception to this is the 
DP PDF: to estimate the systematic uncertainty from 
the mES sideband extrapolation, we use large samples of 
e + e~ — * qq MC data (q — u,d,s,c). We compare the 
distributions of m! and 9' between sidebands at different 
ranges in tties and find the two such sidebands that show 
the maximum discrepancy. We assign as systematic un- 
certainty the effect seen when weighting the continuum 
DP PDF by the ratio of these two data sets. 

The uncertainties associated with Amd and t are esti- 
mated by varying these parameters within the uncertain- 
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FIG. 9: Two-dimensional scans of — 2A log £ as a function of (3 e g and C (top) and the one-dimensional scans as a function of 
/3 e ff (bottom) for the fo(980)K$ (left) and p°(770)Kg (right) isobar components. The value — 2Alog£ is computed including 
systematic uncertainties. On the two-dimensional scans, shaded areas, from the darkest to the lightest, represent the one to five 
standard deviations contours. The statistical (dashed line), and total (solid line) — 2 A log £ are shown on the one- dimensional 
scans, where horizontal dotted lines mark the one and two standard deviation levels. 
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FIG. 10: Two-dimensional scans of -2Alog£ as a function of (S,C), for the f (980)K^ (left) and p°(770)K° s (right) isobar 
components. The value — 2Alog£ is computed including systematic uncertainties. Shaded areas, from the darkest to the 
lightest, represent the one to five standard deviations contours. The • (*) marks the expectation based on the current world 
average from b — > ccs modes [j| (zero point). The dashed circle represents the physical border S 2 + C 2 = 1. 
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A<MK*(892)jt) WW/) 

FIG. 11: Statistical (dashed line) and total (solid line) scans of — 2Alog£ as a function of the relative phases A&(K* (892)-7r) 
(left) and A&((Kir)o) (right). Horizontal dotted lines mark the one and two standard deviation levels. 




FIG. 12: Statistical (dashed line) and total (solid line) scans of — 2Alog£ as a function of the phase differences 
c/>(K* (892)n, {K-kYq'k) (left), 0(/o (980) p° (770)7^) (middle), and </>{p {770)K° s , iT(892)7r) (right). The top (bottom) row 
shows B° (B°) candidates. Horizontal dotted lines mark the one and two standard deviation levels. 




0.005 0.01 0.015 0.02 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.05 0.1 0.15 0.2 



FF(X C0 K°) FF(f x (1300)K°,f 2 (1270)K°) FF(NR) 

FIG. 13: Statistical (dashed line) and total (solid line) scans in terms of — 2Alog£ as a function of the fit fractions of the 
Xco(lP)Kg component (left), the sum of fit fractions of the /2(1270)A'g and fx(^300)Kg components (center), and the flat 
phase space NR component (right). These scans are used to extract the probability of null values of these fit fractions. 
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FF(total) a£°' 

FIG. 14: Statistical (dashed line) and total (solid line) scans of — 2A log C as a function of the total fit fraction (left) and the 
inclusive direct CP-asymmetry A'^p 1 (right). A horizontal dotted line marks the one standard deviation level. 




FF(K*(892)ji) FF((Kn)* o it) 

FIG. 15: Statistical (dashed line) and total (solid line) scans of -2Alog£ as a function of the fit fractions FF(f (980)K%) (top 
left), FF(p°(770)Ks) (top right), FF(A"* ± (892)tt t ) (bottom left), and FFUKiv)^^) (bottom right). A horizontal dotted 
line marks the one standard deviation level. 



ties on the world average [3l| . 

The signal PDFs for the At resolution and tagging 
fractions are determined from fits to a control sample of 
fully reconstructed B decays to exclusive final states with 
charm, and the uncertainties are obtained by varying the 
parameters within the statistical uncertainties. 

Finally, the uncertainties due to particle identification, 
tracking efficiency corrections, Kg reconstruction, and 
the calculation of N BS are 2.0%, 1.6%, 0.9%, and 1.1%, 



respectively. These contribute only to the branching frac- 
tion systematic uncertainties. 

The average fraction of misreconstructed signal events 
(/scf) predicted by the MC simulation has been verified 
with fully reconstructed B — > Dp events [42]. No signif- 
icant differences between data and the simulation were 
found. To estimate a systematic uncertainty from /gcF; 
we vary these fractions, for all tagging categories. Tag- 
ging efficiencies, dilutions, and biases for signal events 
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A cp (K*(892)n) A cp ((Kji)* o ji) 

FIG. 16: Statistical (dashed line) and total (solid line) scans of — 2Alog£ as a function of the direct CP asymmetries 
Acp(K* ± (892) 7r T ) (left) and Acp^K-n)^^) (right). A horizontal dotted line marks the one standard deviation level. 
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FIG. 17: Distributions of At when the B^ ag is a B° (top), B° (middle), and the derived At asymmetry (bottom) for events in 
the J/tp Kg region. The solid line is the total PDF and the points with error bars represent data. 
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TABLE VII: Summary of systematic uncertainties on Q2B parameters. Errors on relative fractions (/3 e ff and phases) are given 
in percent (degrees). 



Parameter 


DP Model 


Lineshape 


Fit Bias 


B Background 


Other 


Total 


C(/„(980)Kg) 

FF(f (980)K° 3 ) 

/3eff(/o(980)Xg) 


0.04 
0.6 
2.1 


0.02 
0.69 
1.9 


<0.01 
0.5 
<0.1 


0.01 
0.07 
0.2 


0.02 
<0.01 
0.3 


0.05 
1.03 
2.9 


C(p°(770)K° s ) 
FF(p° (770) K%) 
0cs(p o (77O)K O c;) 


0.03 
0.23 
1.8 


0.04 
0.31 
2.2 


<0.01 
0.3 
<0.1 


0.06 
0.09 
1.2 


0.06 
0.15 
1.7 


0.10 
0.52 
3.5 


Acp(A-* (892)tt) 

FF(Jf*(892)7r) 

A<E>(7y*(892)tt) 


0.02 
0.8 
8.1 


<0.01 
0.13 
2.8 


<0.01 
0.4 
<0.1 


<0.01 
0.03 
1.4 


<0.01 
0.43 
3.3 


0.02 
1.00 

9.3 


A C p((A:vr) ( » 

FF((^tt) ( » 

A$((^) ( » 


0.02 
0.90 
4.4 


<0.01 
0.39 
2.4 


<0.01 
1.8 
<0.1 


<0.01 
0.12 
1.3 


0.02 
0.33 
3.0 


0.03 
2.08 
6.0 


C(/ 2 (1270)Kg) 
FF(/ 2 (1270)Xg) 


0.07 
0.69 


0.04 
0.16 


<0.01 
0.09 


0.05 
0.02 


0.06 
0.19 


0.11 
0.74 


C7(/x(1300)Jf§) 
FF(/x(1300)Xg) 


0.09 
0.87 


0.03 
0.28 


<0.01 
0.14 


0.01 
0.02 


0.03 
0.17 


0.10 
0.94 


C(jVi?) 
FF(NR) 


0.04 
0.60 


0.01 
0.86 


<0.01 
0.5 


0.01 
0.12 


0.07 
1.62 


0.08 
2.00 


C(XooK° s ) 
FF( Xc0 K° s ) 


0.05 
0.09 


0.02 
0.06 


<0.01 
0.04 


0.01 
<0.01 


0.02 
<0.01 


0.06 
0.11 


A CP 

FFrot 


<0.01 
1.15 


<0.01 
0.55 


<0.01 
2.0 


<0.01 
0.08 


0.01 
0.36 


0.01 
2.40 


^(M980)K s,p°(770)Kl) 
^(p°(770)K° s ,K*(892)7r) 
0(p°(77O)^§,(7f7r) ( » 
0(A'*(892)7r,(A:7r) ( » 


4.4 
12.7 
8.7 
4.7 


2.6 
3.0 
8.5 
0.7 


<0.1 
<0.1 
<0.1 
<0.1 


3.4 
3.6 
3.9 
0.3 


4.3 
7.3 
3.7 
4.6 


7.5 
15.4 
13.3 
6.6 


Signal Yield 


31.7 


5.8 


14.0 


3.3 


23.0 


42.1 
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are varied within their experimental uncertainties. 



VIII. SUMMARY 

We have presented results from a time-dependent 
Dalitz plot analysis of B° — ► Kgir + Tr~ decays obtained 
from a data sample of 383 million T(4f>) — > B5 decays. 
Using an amplitude analysis technique, we measure 15 
pairs of relative magnitudes and phases for the different 
resonances, taking advantage of the interference between 
them in the Dalitz plot. From the measured decay am- 
plitudes, we derive the Q2B parameters of the resonant 
decay modes. Two solutions, with equivalent goodness- 
of-fit, were found. 

Including systematic and Dalitz plot model uncer- 
tainties, the combined confidence interval for the mea- 
sured values of /3 e g in B° decays to / (980)Alj is 18° < 
PcS < 76° at 95% C.L. CP conservation in B° decays 
to fo(980)Kg is excluded at 3.5cr, including systematics. 
For B° decays to p°(770)Ag, the combined confidence 
interval is -9° < /3 cff < 57° at 95% C.L. These results 
are both consistent with the measurements in b — > ccs 
modes. 

In decays to A*+(892)tt-, we find A C p = -0.20 ± 
0.10 ± 0.01 ± 0.02. For the relative phase be- 
tween decay amplitudes of B° — > A* + (892)7r~ and 
B° -> A*-(892)7r+, we exclude the interval -137° < 
A$(A*(892)tt) < -5° at 95% C.L. This last result, com- 
bined with measurements of branching ratios, direct CP 
asymmetries, and relative phases in A* + (892)7r~ and 
A*°(892)7r°, plus a theoretical hypothesis on the contri- 
butions of electroweak penguins to the decay amplitudes, 
can be used to set non-trivial constraints on the CKM 
parameters (~ p, rf ) b y following the methods proposed in 
Refs. [II [II [M El. 
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APPENDIX 

The full (statistical, systematic, and model depen- 
dence) correlation matrices of the isobar parameters for 
solutions I and II are given in Tables I Villi and IIX1 
respectively. The tables are organized in blocks for 
c, c, argc, and argc. Here, the abbreviations fo, p , 
K* , S, /2, fx, NR, and x represent the components 
/o(980)A°, p°(770)A°, A*(892)tt, (Kir) fa / 2 (1270)A°, 
/x(1300)AT<j, nonresonant, and XcoKg, respectively. 
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TABLE VIII: Full correlation matrix for the isobar parameters of solution I. The entries are given in percent. Since the matrix 
is symmetric, all elements above the diagonal are omitted. 





p° 


K* 


S 


|c| 
h 


fx 


NR 


X 


fo 


p° 


K* 


|c 

S 


h 


fx 


NR 


X 




P° 


100.0 
































K* 


51.9 


100.0 






























S 


54.0 


65.0 


100.0 


























1 1 

|c| 


r 

h 


8.4 


2.8 


21.0 


100.0 


























fx 


14.9 


23.2 


32.2 


22.7 


100.0 
























NR 


5.2 


35.0 


24.4 


12.6 


39.3 


100.0 






















X 


6.4 


9.9 


7.8 


2.0 


7.4 


6.1 


100.0 




















fo 


31.3 


30.3 


39.9 


25.2 


36.7 


31.3 


8.0 


100.0 


















n 

p 


20.6 


48.6 


51.2 


8.0 


27.7 


27.5 


5.6 


17.3 


100.0 
















K* 


44.7 


73.5 


56.3 


-4.8 


24.9 


22.0 


9.5 


22.6 


43.4 


100.0 














c< 
O 


59.6 


71.9 


79.7 


21.8 


39.3 


26.9 


11.3 


35.2 


49.4 


57.7 


100.0 










|c| 


h 


2.4 


-10.1 


6.3 


-56.1 


-1.5 


3.9 


-0.3 


10.7 


-6.2 


-21.5 


5.0 


100.0 










fx 


14.5 


34.1 


12.5 


16.1 


-23.0 


12.4 


2.5 


34.5 


7.3 


8.3 


12.9 


-6.2 


100.0 








AT Tt 

N R 


17.8 


57.6 


41.7 


12.7 


10.1 


49.7 


2.4 


40.0 


32.1 


25.0 


31.7 


7.5 


46.2 


100.0 






X 


18.9 


27.0 


30.6 


5.8 


11.8 


9.5 


-84.2 


21.5 


17.8 


24.1 


27.8 


0.8 


8.1 


20.2 


100.0 




P° 


-11.2 


13.3 


4.0 


-16.1 


-2.9 


-2.1 


-0.5 


—0.2 


24.1 


16.3 


3.2 


-3.3 


8.9 


2.1 


4.2 




K* 


25.0 


8.6 


-3.2 


-0.2 


-15.7 


-9.7 


6.3 


-10.4 


-3.9 


5.5 


16.0 


3.8 


6.3 


-6.5 


-3.2 




S 


33.0 


19.6 


3.4 


-4.7 


-17.3 


-16.5 


6.2 


-9.6 


1.0 


18.7 


21.3 


-4.2 


9.6 


-4.2 


1.1 


arg(c) 


h 


12.1 


-0.6 


-9.8 


-2.6 


-23.1 


-27.4 


0.9 


-16.7 


-7.2 


2.2 


1.1 - 


-10.6 


7.2 


-14.1 


-2.6 




fx 


25.0 


10.2 


5.4 


-0.5 


-11.4 


-11.8 


1.0 


-0.8 


2.6 


8.5 


11.8 


-3.8 


15.6 


2.4 


0.4 




NR 


31.6 


17.0 


39.3 


1.0 


-27.1 


-31.7 


-6.7 


11.3 


12.8 


14.5 


19.0 


3.3 


21.5 


19.6 


14.2 




X 


8.6 


1.8 


9.8 


0.6 


-9.9 


-8.9 


-7.9 


2.8 


3.8 


1.3 


4.2 


3.5 


7.3 


8.9 


12.4 




fo 


32.2 


11.7 


18.9 


3.5 


-20.3 


-26.2 


-1.6 


-3.6 


-6.9 


7.3 


18.2 


1.8 


20.3 


-7.1 


4.3 




P° 


14.5 


18.0 


14.6 


-17.3 


-13.4 


-21.0 


-0.7 


-8.7 


14.3 


19.8 


13.4 


1.7 


7.2 


-4.4 


5.4 




K* 


17.1 


7.1 


22.0 


5.2 


-13.5 


-17.3 


-2.1 


5.0 


7.2 


6.5 


13.8 


8.1 


12.8 


29.5 


9.6 


arg(c) 


S 


22.5 


15.9 


25.2 


-3.2 


-16.9 


-21.6 


-0.5 


4.2 


10.6 


17.7 


16.1 


1.7 


14.1 


28.8 


10.8 




h 


15.1 


4.9 


15.5 


-5.0 


-15.5 


-17.9 


-2.1 


10.0 


-2.5 


3.9 


2.9 


11.1 


15.7 


18.6 


7.5 




fx 


8.1 


2.7 


12.3 


-0.6 


16.5 


-20.4 


-0.9 


12.2 


6.1 


3.4 


4.8 


1.4 


-14.6 


4.7 


6.5 




NR 


15.3 


4.1 


14.5 


-3.0 


-22.6 


-20.8 


0.8 


1.7 


8.2 


1.8 


5.2 


2.6 


20.0 


15.1 


3.2 




X 


10.9 


1.1 


12.8 


0.7 


-13.9 


-18.0 


-4.7 


2.1 


3.3 


0.6 


3.9 


5.9 


9.8 


13.4 


8.2 













arg(c) 














arg 


-(c) 










P° 


K* 


S 


h 


fx 


NR 


X 


fo 


P° 


K* 


S 


h 


fx 


NR X 




P° 


100.0 






























K* 


10.4 


100.0 




























S 


18.2 


90.9 


100.0 
























arg(c) 


h 


19.6 


54.1 


61.8 


100.0 
























fx 


25.5 


49.3 


56.9 


58.1 


100.0 






















NR 


24.3 


17.2 


29.9 


31.6 


47.8 


100.0 




















X 


5.0 


6.7 


7.9 


10.2 


17.6 


30.8 


100.0 


















fo 


18.0 


34.3 


42.0 


39.8 


52.9 


55.6 


23.8 


100.0 
















P° 


55.3 


22.2 


32.4 


25.7 


36.6 


42.2 


17.4 


58.8 


100.0 














K* 


4.0 


21.5 


28.0 


23.2 


36.1 


53.9 


31.3 


46.8 


33.5 


100.0 










arg(c) 


S 


9.6 


23.7 


35.1 


27.8 


41.2 


60.7 


33.3 


53.4 


42.7 


90.9 


100.0 










h 


5.5 


6.4 


12.4 


1.5 


29.3 


46.4 


23.5 


44.1 


36.7 


56.7 


60.8 


100.0 








fx 


1.7 


0.0 


5.4 


13.8 


15.5 


36.4 


19.5 


22.2 


22.5 


42.1 


44.8 


39.4 


100.0 






NR 


7.2 


19.2 


27.5 


28.9 


42.3 


55.5 


32.9 


47.3 


37.9 


63.2 


72.5 


48.1 


48.4 


100.0 




X 


4.1 


8.9 


13.3 


15.5 


27.1 


43.3 


35.9 


38.0 


26.9 


55.9 


58.9 


40.0 


33.6 


52.1 100.0 
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TABLE IX: Full correlation matrix for the isobar parameters of solution II. The entries are given in percent. Since the matrix 
is symmetric, all elements above the diagonal are omitted. 









\c\ 










\A 








p° K* 


S 


h 


fx NR 


X 


fo 


p° K* 


s h 


fx NR 


X 



p 

K* 
S 

h 

fx 
NR 

X 



100.0 
46.9 
49.1 

8.7 
16.8 
-8.4 

5.5 



100.0 
68.2 
7.7 
40.3 
30.2 
11.7 



100.0 

25.4 
38.5 
21.2 
9.3 



100.0 
26.6 
9.4 
3.4 



100.0 
49.9 
12.1 



100.0 
9.1 



100.0 





fo 


29.2 


42.1 


50.2 


31.5 


57.9 


34.1 


10.0 


100.0 


















P° 


61.5 


68.1 


40.4 


6.9 


20.6 


6.4 


6.0 


31.6 


100.0 
















K* 


39.8 


75.7 


59.8 


0.3 


33.1 


25.3 


10.9 


33.2 


36.3 


100.0 














S 


50.6 


75.2 


83.2 


25.4 


49.9 


33.4 


13.1 


51.6 


46.0 


61.4 


100.0 












h 


0.8 


-6.1 


9.6 


-53.9 


6.0 


13.3 


0.2 


14.7 


5.3 


-18.5 


10.4 


100.0 










fx 


10.0 


-3.3 


-0.9 


-10.6 


-68.7 


-17.8 


-5.2 


-18.4 


6.3 


-4.0 


-4.9 


2.2 


100.0 








NR 


23.1 


68.8 


44.7 


13.5 


39.3 


34.4 


5.8 


45.6 


58.3 


32.8 


45.4 


14.7 


-13.8 


100.0 






X 


22.3 


33.5 


37.8 


9.8 


19.3 


9.9 


-79.2 


31.3 


20.7 


30.2 


36.1 


3.3 


-2.6 


23.3 


100.0 




P° 


-23.1 


13.7 


5.5 


-11.4 


8.0 


5.2 


0.0 


9.0 


-11.9 


14.5 


6.3 


-0.2 


0.3 


3.8 


6.9 




K* 


30.6 


2.0 


-2.2 


-6.3 


-16.1 


-28.1 


-0.0 


-15.2 


14.3 


-1.4 


6.1 


2.0 


19.4 


-10.3 


0.5 




S 


38.1 


8.9 


1.8 


-10.1 


-17.9 


-39.5 


-0.1 


-15.8 


17.4 


9.4 


7.7 


-8.2 


19.7 


-12.1 


3.7 


arg(c) 


h 


18.1 


-10.0 


-13.7 


-7.4 


-15.4 


-41.3 


-2.4 


-18.6 


1.0 


-6.2 


-10.2 


-12.7 


10.7 


-21.6 


-2.7 




fx 


26.2 


-7.8 


-12.2 


-5.9 


-7.7 


-35.9 


-1.7 


-14.5 


7.8 


-5.7 


-8.8 


-9.9 


12.2 


-15.2 


-3.6 




NR 


32.4 


-0.4 


21.4 


0.5 


-29.5 


-65.2 


-10.4 


-4.2 


12.0 


0.2 


0.4 


-6.4 


21.2 


-8.1 


10.2 




X 


15.4 


-2.2 


0.2 


-1.6 


-9.9 


-18.3 


-4.9 


-5.6 


5.6 


-3.0 


-0.2 


-0.8 


9.2 


-5.6 


4.0 




fo 


30.1 


-8.0 


-2.3 


-0.9 


-13.2 


-43.0 


-2.8 


-16.7 


12.1 


-7.2 


-5.5 


-4.9 


10.4 


-18.7 


-1.6 




P° 


7.6 


11.4 


5.8 


-7.5 


-1.8 


-24.7 


0.6 


-7.5 


4.1 


15.1 


5.5 


-12.6 


1.3 


-7.0 


4.0 




K* 


27.0 


0.8 


7.6 


5.6 


2.8 


-27.8 


0.6 


-2.0 


9.1 


1.5 


7.1 


3.2 


-6.9 


13.9 


4.1 


arg(c) 


S 


32.6 


8.0 


8.4 


-1.1 


0.6 


-31.3 


2.1 


-4.1 


12.6 


12.1 


7.6 


-5.6 


-4.4 


12.2 


4.7 




h 


18.7 


1.7 


6.6 


10.1 


9.8 


-22.9 


0.7 


7.6 


8.6 


3.5 


0.6 


-5.6 


-21.6 


9.3 


4.6 




fx 


21.9 


1.8 


4.4 


9.6 


-0.7 


-30.2 


0.1 


-5.0 


8.1 


2.8 


4.0 


-17.3 


1.0 


-6.6 


-0.2 




NR 


27.7 


-1.9 


-3.0 


3.9 


-0.5 


-30.7 


2.8 


-13.3 


7.8 


-1.5 


-1.2 


-13.8 


-7.2 


-3.7 


-5.0 




X 


19.7 


-5.0 


-0.5 


2.3 


-4.4 


-27.6 


2.7 


-6.1 


6.2 


-4.1 


-2.5 


-1.6 


-0.2 


-0.1 


-2.9 













arg(c) 














arg 












P° 


K* 


S 


h 


fx 


NR 


X 


fo 


P° 


K* 


S 


h 


fx 


NR X 




P° 


100.0 






























K* 


2.9 


100.0 




























S 


7.4 


90.6 


100.0 
























arg(c) 


h 


9.9 


56.6 


65.5 


100.0 
























fx 


5.9 


57.0 


64.4 


69.5 


100.0 






















NR 


10.1 


37.0 


50.3 


44.4 


46.6 


100.0 




















X 


2.6 


39.3 


40.3 


29.1 


31.3 


28.6 


100.0 


















fo 


-0.6 


45.8 


53.5 


47.1 


61.0 


51.9 


27.4 


100.0 
















P° 


41.3 


29.5 


39.2 


31.2 


39.1 


33.0 


16.5 


54.9 


100.0 














K* 


-11.6 


35.2 


39.7 


30.4 


42.7 


30.0 


17.6 


56.0 


32.9 


100.0 










arg(c) 


S 


-8.7 


38.8 


47.7 


36.1 


49.1 


33.7 


19.5 


62.4 


41.1 


91.1 


100.0 










h 


-5.4 


12.2 


17.9 


7.0 


28.5 


27.2 


9.9 


52.8 


42.0 


59.3 


61.6 


100.0 








fx 


-7.0 


23.2 


28.6 


28.0 


34.4 


29.9 


15.4 


34.6 


30.2 


43.3 


47.1 


41.5 


100.0 






NR 


-9.0 


41.4 


47.9 


44.2 


59.5 


30.9 


25.2 


68.9 


48.1 


68.6 


77.6 


54.6 


55.8 


100.0 




X 


-7.3 


29.3 


33.3 


28.8 


38.8 


29.9 


8.8 


47.1 


26.9 


54.7 


58.0 


38.6 


35.8 


54.6 100.0 
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